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ABSTRACT 

Use of a wave rotor as a topping cycle for a gas turbine engine 
can improve specific power and reduce specific fuel 
consumption. Maximum improvement requires the wave rotor 
to be optimized for best performance at the mass flow of the 
engine. The optimization is a trade-off between losses due to 
friction and passage opening time, and rotational effects. An 
experimentally validated, one-dimensional CFD code, which 
includes these effects, has been used to calculate wave rotor 
performance, and find the optimum configuration. The technique 
is described, and results given for wave rotors sized for engines 
with sea level mass flows of 4, 26, and 400 lb/sec. 

INTRODUCTION 

Wave rotors are devices which use unsteady waves to compress 
and expand air in a single device. Basically, a wave rotor 
consists of a rotating bank of passages, inside a housing, with 
ports bringing air in, or taking air out. A four port wave rotor is 
shown in figure 1. Four port wave rotors have been tested by 
General Electric [1], and Rolls -Royce [2]. Other arrangements 
are possible [3]; however, the four port scheme is more attractive 
because it fits easily into a conventional engine. Consequently, 
only the four port scheme will be considered here. With a wave 
rotor installed as a topping cycle in a gas turbine engine, the air 
from the engine compressor would be directed into the wave 
rotor at the input port 1. (see figure 1). The air flows into the 
passages on the rotor, and is compressed by a series of 
compression or shock waves. This air leaves the wave rotor from 
port 2, at a higher pressure than when it entered the wave rotor, 
and passes to the burner. After being heated in the burner, the 
gas returns to the wave rotor through port 3, driving a shock into 
the air in the passages in the process. This gas is trapped on 
board as port 3 closes, at quite high pressure. When the rotor 
comes around to the exit port 4, the gas expands out into the 
lower static pressure in the port. From there the gas goes to the 
engine turbine. In passing through the wave rotor, the air is first 
compressed and then expanded. Thus the wave rotor combines 
in a single device the functions performed by the compressor and 
turbine in a high spool. By using a wave rotor topping cycle, 
combustion temperatures greater than the turbine inlet 
temperature can be used, since the gas leaving the combustor is 
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Figure 1 Schematic of a four port wave rotor 

cooled in the expansion before being sent to the turbine. Also, 
since the rotor is washed alternately by cool inlet air and hot 
combustion gas, it is self-cooled, and attains a steady state 
temperature significantly lower than the combustion temperature. 
By increasing the overall cycle pressure ratio, and allowing 
higher combustion temperatures, the wave rotor topping cycle 
offers a potential route to higher engine efficiency. 

The successful application of a wave rotor topping cycle is 
dependant upon how much performance improvement it can 
achieve over a conventional cycle. Calculations of 
improvements in specific power and efficiency have been 
reported [4]. In making such calculations, it is necessary to know 
the performance of the wave rotor. What is meant by wave rotor 
performance is the pressure gain created by the wave rotor as a 
function of the temperature ratio across it, i.e. P/Pj versus T 4 /T,. 
Note that this pressure gain, P 4 /Pj, is equivalent to the pressure 
at the exit of a turbine divided by the pressure at the entrance to 
the compressor. It should not be confused with the compressor 
pressure ratio, which is commonly used as a measure of 
performance in gas turbine engine cycles. In [4], the wave rotor 
pressure gain was taken to be that measured by Klapproth et 
al.[l]. This is a conservative assumption, in that this 
performance has clearly been demonstrated, but has certainly not 
been optimized. The better the performance of the wave rotor. 


the greater will be the increase in specific power, and reduction 
in specific fuel consumption of the engine/wave-rotor 
combination. The wave rotor performance will depend on the 
losses, which in turn depend on the geometry, and inlet 
conditions, of the wave rotor. What is needed then, is a realistic 
calculation of wave rotor performance, including the major 
losses, which can be used to examine variations in geometry to 
find optimum performance. A suitable, one-dimensional, 
computational fluid dynamics (CFD) code has been written [5,6], 
and was used for this report. The code has been validated against 
experiment [7], and gives good agreement with the experimental 
results of Klapproth et al. [1], and Kentfield [8]. Using this code, 
an optimum wave rotor size has been found for hypothetical 
engines with mass flows typical of small, medium, and large gas 
turbine engines. Finally, the specific power and specific fuel 
consumption of these engines was calculated assuming a wave 
rotor topping cycle employing the optimized wave rotor, for 
comparison with the basic engine. 

ADVANTAGES OF USING A WAVE ROTOR 

The advantages of using a wave rotor topping cycle can be 
calculated with a simple model [4]. The details of the calculation 
are standard, and will not be repeated here; only an outline will 
be given. Because in the topping stage (i.e. the wave rotor - 
burner combination) the gas is heated in a constant pressure 
combustor, thermodynamically it is a Brayton cycle. Thus, it can 
be represented on a temperature-entropy diagram as in figure 2. 
The numbers in figure 2 correspond to the stations indicated in 
figure 1, i.e, 0 is the entrance to the engine, 1 is the exit from the 
compressor, and entry to the wave rotor, 2 is the high pressure 
exit from the wave rotor, and entry to the combustor, 3 is the exit 
from the combustor, and re-entry to the wave rotor, 4 is exit from 
the wave rotor to the turbine, 5 is exit from the turbine, and 6 is 
return to ambient pressure. The basic engine, to which the wave 
rotor is added, is a simple gas turbine engine, so that the work of 
compression is equal to the turbine work ( except for mechanical 
losses ), and the available energy is simply the difference 
between the energy of expansion from a pressure P t , and 
temperature T 4 at turbine inlet to ambient conditions (at pressure 
P 0 ), minus the turbine work. With a wave rotor topping stage 
added, because the stage acts as a pressure gain combustor, the 
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Figure 2 Temperature-entropy diagram for a gas turbine engine, 
and for a gas turbine engine with a wave rotor topping cycle 


Pressure Ratio across Wave Rotor, P4/P1 



Figure 3 Wave rotor performance measured experimentally by 
Klapproth et al. [1], (•), together with a curve-fit calculated with 
the CFD code ( — ) 

pressure P 4 leaving the wave rotor and entering the turbine is 
higher than the pressure P, entering the turbine in the case of the 
basic engine alone. Consequently the expansion between turbine 
inlet and ambient conditions is through a higher pressure ratio, 
and there is more available work. 

It will be assumed that for the basic engine, the value of the 
compression ratio H = P/P 0 is known, as well as the turbine inlet 
temperature, and the compression and turbine efficiencies. With 
these quantities given, the specific power and specific fuel 
consumption of the basic engine can be calculated. When a wave 
rotor topping stage is added, it will be assumed that the turbine 
inlet temperature remains unchanged, but the pressure entering 
the turbine is raised from P 2 to 5 • Thus to calculate the 
performance in this case, it is necesssary to know the pressure 
ratio P 4 /P i generated by the wave rotor for a given temperature 
ratio across it, i.e. T 4 /Tj. 

The wave rotor pressure ratio as a function of the wave rotor 
temperature ratio is not a unique function. It depends on the 
geometry of the wave rotor, which will change with the mass 
flow of the engine. In [4] a fit to the results of Klapproth et al. 
[1] was used. As explained above, use of this fit is conservative; 
the performance should be amenable to improvement, since the 
experimental rotor was neither well timed, nor optimized. 
Moreover, the experimental rotor had a small mass flow; in 
general wave rotor performance improves as the mass flow 
increases. The curve-fit to the experimental results was obtained 
using the one-dimensional CFD code described below, with 
estimates of the geometry of the G.E. rotor, to calculate the wave 
rotor pressure ratio at different values of the wave rotor 
temperature ratio. The fitted curve is the line shown in figure 3. 
Also shown are the experimental points from [1], which are the 
solid circles. Although the fit to the experimental points appears 
to be quite good, this must be regarded as somewhat fortuitous. 
The timing (i.e. when the ports open and close) of the 
experimental rotor was not known, neither was the exact 
geometry, nor the size of any leakage gaps. Consequently an 
exact calculation was not possible. Nevertheless, the fit is 
obviously of the right magnitude, and should provide a good 
estimate of the wave rotor performance for this rotor. 
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Figure 4 Specific power vs. shaft compression ratio for 
engines with and without a wave rotor topping cycle. The 
engines are assumed to have performance as calculated with 
the fit to the data of Klapproth et al. [1] 
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Figure 5 Specific fuel consumption vs. shaft compression ratio 
for engines with, and without a wave rotor topping cycle. The 
wave rotors are assumed to have performance as calculated 
with the fit to the data of Klapproth et al. [1] 

With the wave rotor pressure ratio known from figure 3, the 
specific power, thermal efficiency, and specific fuel consumption 
can be calculated. In [4], calculations of the thermal efficiency 
of an untopped engine, calculated as above, were compared with 
results from the ONX program [9], and found to give excellent 
agreement, giving confidence in this approach. 

With this model, and the curve-fit to the G.E. data, values of 
specific power and specific fuel consumption have been 
calculated for shaft compression ratios (i.e., the compression 
ratio of the mechanical compressor) between 5 and 60, for 
engines with, and without, a wave rotor topping cycle. Turbine 
inlet temperatures of 2400 R, typical of a small engine, 2800 R, 
and 3200 R, appropriate for an advanced engine, were used. Sea 
level ambient conditions were assumed. The results are shown 
in figures 4 and 5. It can be seen that quite large increases in 
specific power are possible by the addition of a wave rotor 
topping cycle, particularly at low values of the shaft compression 
ratio. For example, at a turbine inlet temperature of 2400 R, and 
a shaft compression ratio of 8, addition of a wave rotor topping 
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Figure 6 Specific fuel consumption vs. shaft compression ratio 
with turbine inlet temperature as a parameter 

cycle increases the specific power by 19.2%, and reduces the 
specific fuel consumption by 16.2%. Raising the turbine inlet 
temperature to 3200 R at a compression ratio of 8 reduces the 
gain. At these conditions, the specific power increases 17.4%, 
and the specific fuel consumption decreases by 14.8%. 
Increasing the compression ratio reduces the gain even more 
dramatically. At a shaft compression ratio of 40, the specific 
power increase is 5.8%, and the specific fuel consumption 
decrease is 5.6%. An engine with a turbine inlet temperature of 
2400 R, and a compression ratio of 40 would actually decrease 
its performance by the addition of a wave rotor topping cycle. 
However, neither an engine with a shaft compression ratio of 8, 
and a turbine inlet temperature of 3200 R, nor an engine with a 
shaft compression ratio of 40, and a turbine inlet temperature of 
2400 R are well designed engines; an actual engine is more likely 
to have a compression ratio somewhere between the value giving 
the minimum specific fuel consumption, and the value giving the 
maximum specific power for the turbine inlet temperature of that 
particular engine. The actual value is strong influenced by the 
mission requirements of the engine. 

In figure 6, the specific fuel consumption of gas turbine engines 
without wave rotor topping cycles, as calculated with the simple 
model described above is plotted against compression ratio, with 
turbine inlet temperature as a parameter, for values of 
compression ratio between those for maximum specific power, 
and those for minimum specific fuel consumption. Also plotted 
are points (i.e. specific fuel consumption and compression ratio) 
for actual engines taken from [10]. In addition, there is a line, 
called the engine design line, which is an approximately average 
line through the data, which will be taken as an appropriate line 
relating compression ratio and turbine inlet temperature for an 
"average" engine. Thus for example, an engine with a 
compression ratio of 20 might be expected to have a turbine inlet 
temperature of 2900 R. Using this relationship between engine 
compression ratio and turbine inlet temperature, the specific 
power and specific fuel consumption were calculated for engines 
with, and without, a wave rotor topping cycle, and the results are 
plotted in figure 7. Again, at low values of shaft compression 
ratio, there is a significant gain; for a shaft compression ratio of 
8, the addition of a wave rotor topping cycle increases the 
specific power by 19%, and reduces the specific fuel 
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Figure 7 Specific power and specific fuel consumption for 
engines with, and without a wave rotor topping cycle. At any 
compression ratio, the basic engines are assumed to have 
turbine inlet temperatures corresponding to the design line of 
figure 6 


consumption by 16.2%. At higher values of shaft compression 
ratio, the advantage is lessened; for a shaft compression ratio of 
40, addition of a wave rotor topping cycle increases specific 
power by 5.8%, and reduces specific fuel consumption by 5.6%. 

It is worth re-emphasizing that the above results were derived 
using the wave rotor performance found experimentally [1]. This 
performance has been demonstrated, and therefore it should be 
possible to obtain the calculated increases in specific power, and 
decreases in specific fuel consumption. However, the 
experimental rotor used in [1] was not optimized even for the 
experimental mass flow, and is certainly not optimum for any 
engine. It will be shown below that wave rotor performance can 
be optimized for a given mass flow, and the resulting 
performance will be better than that demonstrated by Klapproth 
et al.[l]. Consequently the increase in specific power, and 
reduction in specific fuel consumption that can be obtained by 
the use of an optimized wave rotor topping cycle will also be 
greater. In order to perform this optimization, a model that can 
calculate wave rotor performance is required. Such a model is 
described in the next section. 

CFD WAVE ROTOR MODEL 

The computational wave rotor model is a one dimensional code 
which follows a single wave rotor passage as it rotates past the 
various ports of the cycle. The working fluid is assumed to be a 
calorically perfect gas. The details of this code have been 
described in previous publications [5-7] and will not be 
presented here. Although several loss mechanisms are accounted 
for, those of interest here result from the effects of viscosity, 
gradual passage opening (so called finite opening time), and 
mixing of non-uniform velocity profiles in the exhaust port 4. 
For reference, viscous effects are accounted for in the model with 
a source term in the governing momentum equation. The 
strength of the source term is determined by local velocity, 
density and a constant friction factor which depends on the 
passage geometry. It is defined as 
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Since the viscosity does not change much with temperature it is 
assumed constant at a value which is the average of ports 1 and 
4. The hydraulic diameter is defined as 
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The passage opening time is defined as 
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Increases in either finite opening time or friction factor result in 
decreased wave rotor rotor performance. Furthermore, 
examination of EQS (1), (2) and (3) reveals that, aside from an 
increase in rotor angular velocity, reducing the opening time 
increases the friction factor. For the purposes of this study the 
code was modified for use as a design rather than an analysis 
tool. In this capacity the prescribed inputs to the code are the 
friction factor, the ratio of specific heats, the passage opening 
time, the ratio of exhaust port to inlet port stagnation 
temperature, T/Tj, and the so called expansion ratio, €, which 
is the ratio of the static pressure in the exhaust port, p 4 , to the 
assumed uniform pressure in the passage just before the exhaust 
port opens. The code then calculates the inlet stagnation 
pressures and temperatures, outlet static pressures, and port 
locations which insure that the various waves are properly timed 
and that, upon completion of the cycle, the total mass flow from 
outlet ports matches the total flow to the inlet ports. It is noted 
that the code operates in the rotor fame of reference and that as 
such, all stagnation properties discussed in this section of the 
paper are to be understood as rotor relative. The balancing of 
mass flow also insures that the passage completes the cycle in 
nearly the same state at which it began so that the cycle may start 
again. The calculation process is an iterative one with several 
sub-processes. These are shown in flowchart form in Appendix 
A and described briefly below. A complete calculated cycle, 
along with the port nomenclature, is shown in figure 8. This plot 
is a passage pressure contour with time (or circumferential 
position) and distance along the passage as the vertical and 
horizontal axes respectively. It is noted that although the design 
code can calculate the appropriate port timing and relative 
stagnation conditions for all four ports in the cycle, only the inlet 
and exhaust portion of the cycle are used in this study and 
discussed in the paper. The iteration begins at t = 0 with both 
ends of the passage closed and the gas in a state of uniform 
pressure and density and velocity (u = 0). The code then 
commences time integration with the opening of the exhaust port 
4 (figure Al). Integration continues until the expansion fan, 
which was initiated by the exhaust port opening, reflects f om the 
inlet side endwall and returns to the exhaust port. When the 
expansion wave reaches the exhaust port, the outflow velocity 
begins to drop. Eventually, it drops below zero and inflow 
occurs. At this point the first sub-iteration begins (figure A2). 
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Figure 8 Pressure contours for a properly timed four port cycle 


temperature ratio T 4 /T, of the engine. It will be assumed that 
leakage is negligible. There are two reasons for this assumption: 
(1) recent experiments at NASA Lewis Research Center [10] 
have indicated that the effect of leakage on performance is 
negligible for values of the leakage parameter G less than 0.01, 
and (2) work is in progress to develop seals to prevent, or 
significantly reduce, leakage. It is relatively easy to design for 
G < 0.01. 

It will be assumed that m , 13 s , and T 4 /Tj are known for the 
engine. For most of the engines of interest, T/T, has a value 
approximately equal to 2. The wave rotor expansion ratio giving 
the best performance depends on T/Tj. For X /TJ close to 2, 
calculations have shown that the optimum value of the expansion 
ratio is e=0.4. This value was assumed in the calculations 
following. Equating the engine mass flow to that of the wave 
rotor gives; 


This iteration locates the exhaust port closing time such that 
when the closing process is complete, the port velocity has just 
reached zero. This closing time is then used for the remainder of 
the major iteration. It is noted that fixing this time is not entirely 
correct since the opening of the inlet port, to be described below, 
alters the behavior of the expansion fan and thus changes the 
time at which the exhaust velocity goes to zero; however, the 
alteration is generally very small, as is the error introduced. 
With the exhaust port timing set and the exhaust process 
complete, the exhaust mass flow and stagnation properties can be 
calculated. From the calculated value of T 4 and the specified 
value of T 4 /T, the inlet temperature T, may be obtained. An 
initial guess at the inlet stagnation pressure Pj is then made (i.e. 
Pj = P 4 ), and the iteration restarts from t= 0. Integration 
continues until the pressure at the inlet end of the passage, 
adjacent to the wall, matches the inlet stagnation pressure P,. At 
this point the inlet port is opened and the iteration continues until 
the flow in the inlet port reverses. Here, the second sub-iteration 
begins (figure A3). This process insures that the inlet port 
completes the closing process just as the shock, which has 
coalesced from the exhaust port 4, reaches the inlet end of the 
passage. The arrival of the shock is detected by either a sudden 
reversal of the inlet flow or a sudden change in passage pressure 
adjacent to the inlet. With the inlet closing complete, the inlet 
mass flow is known. This is compared to the exhaust mass flow. 
If they match, the iteration is complete, the correct port timing is 
known, and the ratio P^j may be calculated. If not, a new guess 
is made for Pj and the calculation is restarted. Once two guesses 
have been made for P„ it may subsequently be updated using 
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where e is the difference between inlet and exhaust mass flow. 


OPTIMIZATION 

Optimization means finding the wave rotor geometry (i.e. L, R, 
b and h) plus rotational speed, which gives the best performance 
when used as the topping cycle for a particular engine. The wave 
rotor is thus required to handle the mass flow of the engine, at the 


rh = pUA 


(5) 


where p and U are calculated at the wave rotor inlet port In fact, 
p and U can not be calculated at the inlet port, because they are 
set by the conditions in the wave rotor, and are found by solving 
the one-dimensional CFD code. But values of p and U are 
required to calculate the wave rotor geometry, which must be 
known to evaluate input parameters for the code. This dilemma 
is solved either by estimating the values for p and U, or using 
values from a previous solution, finding a solution with these 
values, and then using the new values from the solution to 
calculate new input. This is repeated until the input values are 
equal to the output values. The convergence is quite rapid; 
usually only three or four iterations are needed. Each of the n 
inlet ports has an area given by 


A n = 4> in Rh = £ (6) 

The value of <t> ta comes from the cycle diagram (figure 8) for the 
particular value of T 4 /T, 
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Thus, if a value of R is chosen, the value of h is determined via 

/ T 

h - m cycle 

2itRpU^ x in 

With h fixed, a value of b is assumed, which fixes the hydraulic 
diameter from EQ (2). Given D h , and assuming a value for L, the 
friction factor (EQ (1», and dimensionless passage opening time 
(EQ (3)) can be calculated. One more quantity, the relative 
temperature ratio, (T 4 /T 1 ) rcl is required for input into the code. 
The one-dimensional CFD code solves the unsteady flow within 
the rotor in the relative frame of the rotor. In this frame, the 
temperature ratio across the rotor will be different from the 
temperature ratio in the absolute frame. From the definition of 
relative stagnation temperature [12], it follows that 



5 


T = T,* 1 + — 

2C o 

Pi 


„ _ Trc i _ W 2 

4 4 ^7 

P4 


(9) 

( 10 ) 


The value of W depends on the rotor geometry, increasing with 
increasing R, but decreasing with increasing L 
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optimized with respect to b , at each value of R, is plotted in 
figure 9 for n= 2, and a mass flow of 4 lb/sec. The relative 
performance is also shown in figure 9. At small values of R, the 
value of W is low, and the absolute performance approaches the 
relative performance, whereas at large values of R, there is 
significant difference between the two. Large values of L reduce 
G), and therefore W, but also increase the friction factor, reducing 
performance. From this plot, values of L and R, and thus b and 
h can be selected which will give optimum performance. For the 
case shown, these values are L = 5.25”, R=2.5", /i=0.83", and 
6=0.28". 


APPLICATION TO DIFFERENT ENGINES 


With W known, and knowing the value of TJTi required, it is 
possible to find the value of (T 4 /T 1 ) rcl . The value of J* 1 
is used in the CFD code to solve for (P ^) 1 * 1 , which is then 
corrected to the absolute value via the isentropic relationship 



( 12 ) 


Large values of W are undesireable, since then the absolute 
performance drops significantly. 

The optimization procedure is as follows: for given values of 
mass flow, and absolute temperature ratio T/Tj, corresponding 
to the engine selected, values of n , L and R are picked, thereby 
determining h. The absolute performance is calculated for 
various values of b. At small values of 6 , the opening time is 
small, but the friction factor is high, reducing performance. At 
large values of 6 , friction is low, but the large opening time 
reduces performance. There is an optimum value of b in 
between, which is selected. This calculation is repeated at 
different values of R, for the same value of L. Then a new value 
of L is chosen, and the process repeated. The performance 
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Figure 9 Relative and absolute wave rotor pressure ratio vs. 
length and radius for wave rotors sized for a mass flow of 4 
lbm/s. The passage width has been optimized at each radius and 
length 


Using the methodology outlined above, the dimensions of 
optimal wave rotors can be found for gas turbine engines of 
various sizes. Such calculations have been made for hypothetical 
engines of small (4 Ib/sec), medium (26 lb/sec), and large (400 
lb/sec) engines. For these engines, typical values of compression 
ratio and turbine inlet temperature were selected, as well as 
estimates of the engine dimensions. The values are listed in table 
1 . Also listed are the dimensions and performance of the optimal 
wave rotor for each engine. It will be seen that the wave rotor 
can be relatively small compared with the engine itself. 

In one sense, there is not in fact an optimum wave rotor at all! 
This is because there is an optimum for each value of n, the 
number of cycles per revolution. Although the highest 
performance comes from n = 1 , the change in performance as n 
is increased is very small. On the other hand, the length of the 
optimum rotor scales almost inversely with n. This is illustrated 


Table 1 Optimized wave rotor specifications 


Engine 

Small 

Med. 

Large 

Core Mass Flow (lbm/s) 

4 

25 

400 

Shaft Compression Ratio 

7 

8 

40 

Turbine Inlet Temp. (R) 

2400 

2400 

3200 

Diameter (in.) 

24 

24 

50 

Length (in.) 

40 

48 

190 

Wave Rotor 

Diameter (in.) 

6 

11 

24 

Length (in.) 

6 

12 

12 

Passage Height (in.) 

0.673 

2.50 

4.51 

Passage Width (in.) 

0.30 

0.48 

0.38 

Rotor Speed (r.p.m.) 

15600 

7700 

4200 

T/T, 

2.3 

2.3 

1.9 

P 4 /P 1 

1.232 

1.275 

1.204 

Performance 

untopped Spec. Pwr. (HP/lbm/s) 

121 

151 

276 

topped Spec. Pwr. (HP/lbm/s) 

149 

183 

300 

% Increase from Wave Rotor 

23 

21 

8.7 

untopped SFC (lbm/HP/hr) 

0.650 

0.529 

0.338 

topped SFC (lbm/HP/hr) 

0.528 

0.436 

0.311 

% Decrease from Wave Rotor 

18.8 

17.5 

8.0 
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Figure 10 Absolute pressure ratio vs. length and radius for wave 
rotors sized for a mass flow of 400 lbm/s. Curves are given for 
two wave rotor cycles per revolution, (n= 2), and for four wave 
rotor cycles per revolution, (n=4) 


in figure 10, showing performance curves for the large engine for 
n=2, and n=A. There is almost no change in performance, but the 
length of the optimal n=A rotor is half that of the n= 2 rotor. 
From a size point of view, the n = 4 rotor seems preferable. 
However, the higher value of n leads to much more ducting, 
which could be a severe complication. In table 1, it was 
accn nwl that 71=2 would be appropriate for the small and medium 
engines, but that n = 4 would be more appropriate for the large 
engine. Obviously, for a real engine, the choice of n will be 
made by the engine designer based on engineering 
considerations. Another complication is indicated in figure 9, in 
which approximate values of the wave rotor hub-tip ratio are 
given. Performance optimization for the small engine drives the 
choice of wave rotor radius to small values, which results in low 
values of the hub-tip ratio. The one-dimensional code does not 
presently account for hub-tip ratio effects. Consequently, the 
lower limit of hub-tip ratio for a wave rotor is not known. 
Modification of the code is required to include this effect. 
Turbomachinery typically has a lower limit of hub-tip ratio of 
0.7. It seems reasonable to use a similar criterion for wave rotors 
also, and this may limit the value of radius used. 

The calculated performance of the three optimal wave rotors 
has been plotted in figure 11. It is seen that each wave rotor is 
predicted to give performance better than that of the G.E. wave 
rotor. This is reflected in the wave rotor enhanced performance 
of the three hypothetical engines, as listed in table 1. All three 
engines will give better performance than that indicated by the 
calculation based on the GE wave rotor performance, increasing 
the specific power another 2.3% for the small engine, 3.1% for 
the medium engine, and 2.7% for the large engine, and reducing 
the specific fuel consumption a further 1.6% for the small engine, 
2.6% for the medium engine, and 2.4% for the large engine. The 
increase in specific power, and reduction in specific fuel 
consumption that can be expected for engines corresponding to 
the engine design line of figure 6, when enhanced by wave rotors 
optimized for the small and large engines is indicated in figure 
12 . 
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Figure 1 1 Performance of wave rotors optimized for the small, 
medium, and large engines together with the performance of the 
G.E. wave rotor measured by Klapproth et al. [1] 



Figure 12 The improvement in gas turbine engine 
performance that can be obtained with a wave rotor topping 
cycle employing wave rotors of differing performance, for 
engines corresponding to the engine design line 

SUMMARY AND CONCLUSIONS 

The advantages of using a wave rotor topping cycle on a gas 
turbine engine have been reviewed. A one-dimensional, 
experimentally validated, CFD model for calculation of wave 
rotor performance has been described briefly. This model was 
used to find optimally sized wave rotors for three hypothetical 
engines. It is concluded that use of a wave rotor topping cycle on 
a gas turbine engine can: (1) Based on demonstrated wave rotor 
performance, give increases in specific power of 19.2% at a 
compression ratio of 8, turbine inlet temperature of 2400 R, and 
5.8% at a compression ratio of 40, turbine inlet temperature of 
3200 R. (2) Based on demonstrated wave rotor performance, 
give reductions in specific fuel consumption of 16.2% at a 
compression ratio of 8, turbine inlet temperature of 2400 R, and 
5.6% at a compression ratio of 40, turbine inlet temperature of 
3200 R. (3) With optimally sized wave rotors, give increases in 
specific power of 21% at a compression ratio of 8, turbine inlet 
temperature of 2400 R, and 8.7% at a compression ratio of 40, 
turbine inlet temperature of 3200 R. (4) With optimally sized 
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wave rotors, give reductions in specific fuel consumption of 
17.5% at a shaft compression ratio of 8, turbine inlet temperature 
of 2400 R, and 8% at a shaft compression ratio of 40, turbine 
inlet temperature of 3200 R. In addition, optimal wave rotors are 
reasonably sized for this application. 

NOMENCLATURE 

A = total wave rotor inlet area from the compressor 
A^ = area of a single wave rotor inlet (= A/n) 
a* = reference speed of sound (in passage just prior to opening 
to exhaust) 

b = wave rotor passage width 

c pl = specific heat of air before entering the combustor 
Cp 4 = specific heat of gas after leaving the combustor 
D h = passage hydraulic diameter 

e = difference between inlet and exhaust mass flow rate used in 
EQ (4) 

G = leakage parameter (= 2 6/h) 
h = wave rotor passage height 
L = wave rotor passage length 
m = engine mass flow 
n = number of cycles on the wave rotor 
Pi = stagnation pressure at port i 
Pi = static pressure at port i 
R = rotor radius at mid-passage height 
Tj = stagnation temperature at port i 
t = time 

U = axial flow velocity at input to rotor 

u = velocity inside a passage 

W = circumferential velocity of rotor 

y= specific heat ratio for combustion gases 

8 = spacing between rotor and end-wall 

e= expansion ratio, ratio of p 4 to pressure in passage before 

opening to the exhaust port 4 

p= viscosity 

p= density at input to rotor 

p*= reference density, i.e. density in passage just prior to 
opening to exhaust port 4 

II s =shaft compression ratio, i.e. P,/P 0 , which is the compression 
ratio of the mechanical compressor 
t= non-dimensional time = t a*/L 
Tcydc= non-dimensional time for one cycle 
t in = non-dimensional time inlet port is open 
non-dimensional passage opening time 
4> in = angle subtended by inlet port 
( 0 = rotor angular velocity (radians/sec) 
subscripts 1, 2, 3, 4= wave rotor ports indicated in figure 1 
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APPENDIX A: CFD ITERATION FLOWCHARTS 


u0<0?]^> — — — | Continue calculation | 


YES 




Figure A1 Four port iteration flow chart 



Figure A2 Exhaust port closing sub-iteration flow chart 
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